%kfold RBF 
clc
clear all
load Gammapre_Noneinput
load output_YBOCS_None
load idx_GammaNone3
input=Gammapre_Noneinput(:,[8,3,4]);
input=input';
output=output_YBOCS_None;

K = 12;
clear net
clear Y perf
CrossValIndices = crossvalind('Kfold', 12, K);
for i = 1:12
    display(['Cross validation, folds ' num2str(i)])
    IndicesI = CrossValIndices==i;
    TempInd = CrossValIndices;
    TempInd(IndicesI) = [];
    xTraining = input(:,CrossValIndices~=i);
    tTrain = output(:,CrossValIndices~=i);
      xTest = input(:,CrossValIndices ==i);
      tTest = output(:,CrossValIndices ==i);

 net = newrb(xTraining,tTrain,1.5,0.01,25,12);
 
Y(i,:) = sim(net,xTest);
tTest1(i,:)=tTest;
y1=round(Y);
perf(i,:) = mse(net,tTest1(i,:),y1(i,:))
%  clear net
end

rmse=sqrt(perf');
display(['Y Model=  ' num2str(y1')]) 
display(['Target=   ' num2str(tTest1')]) 
disp(['Expected Error= ',num2str(max (output)*0.1)])
display(['MSE= ' num2str(perf')]) 
display(['Total MSE= ' num2str((sum(perf)/12))]) %mean(perf)
display(['RMSE= ' num2str(rmse)]) %mean(perf)
display(['mean RMSE= ' num2str(mean(rmse))]) %mean(perf)